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Abstract 

We investigate the Dirichlet-scalar equivalent of Casimir-Polder forces between 
an atom and a surface with arbitrary uniaxial corrugations. The complexity of the 
problem can be reduced to a one-dimensional Green's function equation along the 
corrugation which can be solved numerically. Our technique is fully nonperturbative 
in the height profile of the corrugation. We present explicit results for experimentally 
relevant sinusoidal and sawtooth corrugations. Parameterizing the deviations from 
the planar limit in terms of an anomalous dimension which measures the power-law 
deviation from the planar case, we observe up to order-one anomalous dimensions at 
small and intermediate scales and a universal regime at larger distances. This large- 
distance universality can be understood from the fact that the relevant fluctuations 
average over corrugation structures smaller than the atom-wall distance. 

1 Introduction 

Casimir forces [Tj between mesoscopic or macroscopic objects as well as Casimir-Polder 
forces [2] between an atom and a surface can be attributed to a reordering of fluctuations 
in the quantum vacuum. Recent years have witnessed considerable progress in measuring 
these forces j3[ HI [5], [6] , paving the way for future application in micro- and nanomechanical 
engineering and single-atom manipulation. 

Particularly for such applications, standard caicuiational techniques for simple flat sur- 
faces are insufficient and a profound understanding of the influence of geometry on these 
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quantum forces is required. In fact, since fluctuations occur on all momentum or length 
scales, quantum forces are strongly affected by the global properties of a given system. 
From a technical perspective, global properties such as geometry or curvature dependencies 
generally require a full understanding of the fluctuation spectrum in a given configuration 
and cannot be dealt with by perturbative expansions with respect to a small geometry 
parameter. 

Therefore, a variety of new field-theoretical methods for understanding and comput- 
ing fluctuation phenomena have been developed in the past few years, such that early 
phenomenological recipes, such as the proximity force approximation (PFA) [7j have been 
overcome by now. Apart from exact results in certain asymptotic limits [HI E] and more 
controllable approximation techniques JTUJ [11], field-theoretical worldline methods have 
lead to efficient algorithms for Casimir energies [T2l [TBI [Ti] . In addition, approaches based 
on scattering theory have proved most successful for finding new exact solutions and effi- 
cient computation schemes [TJ)J [TBI El [321 ESI [201 [2H [22] , see [23] for a recent review. New 
results with direct mode summation have been obtained in [21]; numerical tools based on 
brute-force discretization have been used in [25] . 

Many of these approaches can be directly linked with a constrained-functional-integral 
formulation, as first introduced in [26J for the parallel-plate case and further developed 
for corrugated surfaces in [27] and lastly extended to general dispersive forces between 
structured media [2H]. For Casimir forces involving corrugated surfaces, results based on 
a perturbative expansion in the height profile have been obtained recently [29|. 

Whereas most of these approaches have been mainly applied to Casimir forces between 
extended objects, the Casimir- Polder force between an atom and a surface of general shape 
has not been so widely studied. An estimate of the influence of surface roughness has been 
given in |30j, where an additivity approximation has been perturbatively expanded in 
the roughness amplitude. In [21], it was pointed out that large geometry corrections to 
the Casimir-Polder force should be observable for atoms in front of a corrugated surface, 
potentially visible in experiments with cold atoms near a surface. These studies were 
based on a perturbative analysis in the height profile. If the height-profile parameters 
are of the same order as the atom-wall distance, however, geometry plays a much more 
dominant role and perturbative theory is expected to be no longer appropriate. Whereas 
cold atom gas experiments become increasingly difficult in this regime, Casimir-Polder 
force measurements can be well controlled in this regime. In particular the atomic beam 
spin echo technique introduced in [321 E3j E] has recently demonstrated high resolution 
access to some of the fundamental issues concerning Casimir interactions. Experimental 
results on geometrical effects will be presented elsewhere shortly [33] . 

This work is devoted to a nonperturbative study of Casimir-Polder forces near a surface 
with uniaxial corrugation. Our method is based on the constrained-functional-integral ap- 
proach which facilitates to map the core part of the problem onto a one-dimensional Green's 
function problem along the direction of nontrivial curvature of the surface. As this Green's 
function problem involves singular kernels, we identify an appropriate representation which 
is accessible to stable and efficient numerical tools. For simplicity, we consider a fluctuat- 
ing scalar field obeying Dirichlet boundary conditions instead of the full electromagnetic 
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field. Hence, our results should not be directly applied to realistic atom-wall studies; they 
are applicable to Casimir configurations in ultracold-gas systems as suggested in [35J. Our 
method is not restricted to Dirichlet scalars and can straightforwardly be generalized to 
the electromagnetic case. 

The paper is organized as follows: In Sect. El we review the treatment of Casimir forces 
within scalar QFT and give the general Casimir energy for two Dirichlet surfaces. Within 
this formulation, we then evaluate the Casimir energy in Sect. [3] for a sphere and a planar 
surface in the limit ^ C 1 to first order explicitly. This limit defines the scalar analogue 
of the Casimir-Polder force. In Sect. HI we extend our technique to the case of uniaxially 
corrugated surfaces. In Sect. we present numerical results for the scalar Casimir-Polder 
potential for a sinusoidally shaped surface as well as for a saw tooth profile, both of which 
have been used in experiments [31]. Conclusions are given in Sect. [71 



2 Scalar field theory with boundaries 

The Casimir energy of a system is given by the shift of the ground state energy caused by 
the presence of boundaries imposing constraints on the fluctuating field. For calculating 
Casimir forces between disconnected bodies, only the Casimir interaction energy is rele- 
vant. The latter corresponds to that part of the ground state shift which depends on the 
relative position and orientation of all bodies. The Casimir self-energies of single bodies 
are irrelevant for the Casimir force. 

Here, we investigate the Casimir interaction energy induced by fluctuations of a scalar 
field 4>. We follow the constrained-functional-integral approach (26], which is introduced 
in the following in a brief and simplified manner, see [261 [2T| I3"9j ED] for a more detailed 
discussion and generalizations. We start from the associated Euclidean generating func- 
tional, 



Z = j P0exp (-±: J d 4 x(d<f>(x)f) 



The Casimir interaction energy of a system bounded by two surfaces, whose relative posi- 
tion is specified (e.g., by a (mean) separation H), is then given by 

E(H) = -^ln^L, (2) 

where Zb.c. stands for the generating functional of the fluctuating field obeying the system's 
boundary conditions, whereas Z^ represents the case of infinite separation between the 
objects, i.e., H — > oo. In this way, irrelevant Casimir self-energies are subtracted. The 
length in Euclidean time direction is denoted by T-&. 

The boundary conditions for the fields are implemented by insertion of a 5 functional 
constraint into Z. For the case of Dirichlet boundary conditions, the corresponding 5 
functional can be represented by product of 5 functions 6(<f>(x a )) for all 4-vectors x a 
pointing onto a surface S a ; here, a labels multiple disjoint surfaces. Hence, Z-q.c. for 
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Dirichlet boundaries in the case of two plates is given by 

Z D = J V<f> f[ H 5 (4(x a )) exp (-1; J d 4 x (<90(x)) 2 ) . (3) 

In order to evaluate the integral over the fields <ft, a Fourier representation is used for the 5 
functional with the help of auxiliary fields that have support only on the surfaces S a \^\ First 
performing the GauBian integral over leaves us - apart from an irrelevant factor - with 
another GauBian integral for the auxiliary fields which can also be carried out, yielding for 
the Casimir interaction energy, 

ftr 1 

E(H) = --Trln(M^M), (4) 
Je * 

where M. denotes a matrix whose entries are the propagators of the fluctuating scalar, i.e., 
the Green's function of the Laplacian in this case, 

M a p(x,x') = —- — (5) 

4tH (x a -Xp) 1 

with four- vectors x and x' pointing onto the surfaces S a and Sp, respectively. At the same 
time, Ai is the inverse of the auxiliary-field propagator on the surfaces. Complications of 
Eq. (j3J) are hidden in the fact that the trace has to be evaluated over the spacetime hyper- 
surface of the boundaries only, as specified by the support of the auxiliary fields. As the 
surfaces can be curved, this requires appropriate metric factors for the spatial integration 
measures (see below). The normalizing operator Aioc corresponds to the propagator at 
infinite separation of the surfaces and takes care of the self-energy subtraction. For the 
case of two surfaces, M.^M. can be written as 

where the matrix components distinguish between spacetime arguments "living" on the two 
different surfaces, i.e., a, (3 = 1, 2. In the matrix product, an integration over the connecting 
intermediate spacetime points is implicitly understood. Due to the normalization with 
Aioo, only the off-diagonal elements are non-trivial, as they contain information about the 
propagation of fluctuations between the two different surfaces. 

In order to compute the Casimir energy from Eq.(j4]) explicitly, we introduce the off- 
diagonal matrix AAi, 

1 Drawing the analogy to the electromagnetic case, the auxiliary fields can be thought of as charged 
sources which enforce the boundary conditions by means of their coupling to the fluctuating field [20 . 
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such that a series expansion of the logarithm in Eq.(j3J) yields 

n=l 

Due to its off-diagonal structure, only even powers of AAi contribute to the Casimir 
energy. Performing the discrete trace over the matrix entries then yields by the cyclicity 
of the trace 

he 1 00 1 

E{H) = £ -Tr ((MTlMnM^M 21 ) n ) . (9) 

n=l 

Thus, the Casimir energy can be understood as a sum over all generic " chains" of correlators 
between the surfaces. Starting from a generic point on surface Si, it is summed over all 
possible auxiliary-field or source correlators on the surface itself, followed by all possible <fr 
propagations to the second surface S2, then generic source propagations on S2 itself and all 
thinkable ways back to S\ again. As a last step, these Casimir contributions are summed 
over all possible starting points. Within the language of scattering theory, is related 
to the T matrix associated with scattering off the surface a [20J . 

Equation flU]) is a useful starting point for explicit computations for a given configura- 
tion. For instance for two parallel Dirichlet planes, the trace turns out to be proportional 
to l/(nH) 3 , yielding the scalar analogue to Casimir's celebrated result upon summation 
over n. The present formalism can straightforwardly be generalized to Neumann or elec- 
tromagnetic boundary conditions by insertion of the corresponding 5 functional constraint. 
The formalism has also been worked out for dielectric materials and fluctuations in media 
from which Lifshitz theory [36J can immediately be derived [28]. Arbitrarily curved sur- 
faces require the knowledge of M-^a f° r both surfaces, which is equivalent to working out 
the T matrix (or S matrix) of the corresponding scattering problem. Moreover, an efficient 
means to carry out the remaining summation and integrations is needed as well. 



3 Scalar Casimir- Polder potential for a planar surface 

The Casimir-Polder force between an atom and a plane wall at a distance H is easily derived 
from Lifshitz theory for two parallel plates. Consider one of the plates as consisting of a 
dilute dielectric with a dielectric permittivity e = 1 + AiraN + 0(N 2 ). Here, a denotes the 
polarizability of the atoms in the surface and N is the number of atoms per unit area in 
the surface. Expanding the Casimir force from Lifshitz theory to first order in iV yields the 
Casimir-Polder force between one of the atoms in the dielectric and the opposite wall as 
the prefactor at order N, see, e.g., [57] . This approach, however, cannot simply be applied 
to the case of an atom near a corrugated surface, since the prefactor at order N is an 
average over all possible atom positions at a given mean distance H above the corrugated 
surface. Hence, the important information about the dependence of the force on lateral 
coordinates is not available in this manner. 
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This information can be extracted following a different strategy. As the Casimir-Polder 
force arises from the atom's polarizability, i.e., the dipole transition induced by the fluc- 
tuating field, we can represent the atom by a compact surface with suitable polarization 
properties. The simplest case is, of course, a small sphere of radius r. In the limit r ^ H, 
the lowest non- vanishing inducible multipole moment, i.e., the dipole polarizability in this 
case, dominates the Casimir interaction between the small sphere and the corrugated sur- 
face. Hence, identifying this dipole polarizability of the sphere with that of an atom yields 
the Casimir-Polder law in this limit. 

In this section, we proceed in exactly this fashion for the Dirichlet scalar case, i.e., we 
define the Casimir-Polder law in this the Casimir force between a small sphere of 

radius r and a corrugated wall at (mean) distance H in the limit of r <C H . Here the 
lowest multipole contribution of a small sphere is that of a monopole excitation, hence the 
resulting power law [13j[T5] with distance H will be different to that of the electromagnetic 
case , but it is straight forward to generalize this approach and computational strategy to 
the electromagnetic case. 

As an exercise, we start with the planar case, identifying the planar surface with S\ 
parametrized by the spatial coordinate 23 = 0. The center of the sphere (surface S2) of 
radius r is located at X\ — x 2 — 0, x 3 — H. The singularity structure of the involved 
propagators is not entirely trivial and special care has to be taken to choose the proper 
order of limits and integrations. Of the four correlators occurring in Eq.^HJ), we begin with 
the contraction of the inverse propagator on the plane with the propagator between the 
surfaces. 

Specifying Eq.Q for the plane Si results in = 1/[4ti 2 {x — x') 2 ], where x = 

{x ,xi,x 2 ) = (xq,x\\), i.e., the x 3 coordinate perpendicular to the plane drops out. By 
Fourier transforming the propagator to postion space, we obtain the functional inverse, 



M n 1 =2^-V 2 , (10) 

where V = {do, 81,82)- This operator can immediately be contracted with the correlator 
between the plane and the sphere yielding 

M^{x : ,d')Mx2{^- 1 x') = - -——^ EE AM12&X'). (11) 

Since the setup is translationally invariant in time direction, it is expedient to Fourier 
transform the propagators to frequency space with respect to the time coordinate. The 
singularity structure can be well controlled using a propertime integral representation (with 
propertime S). In addition, we parametrize the surface of the sphere S2 by spherical 
coordinates. The combined propagator then reads 

AM 12 {C,xfn) = \{H + rcos6) [ dSVSexp (-7-) 

x exp (— [(:ei — r cos sin 6) 2 + {x 2 — r sin (ft sin 9) 2 + {H + r cos #) 2 ] S) , 

(12) 
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where ( denotes the imaginary frequency, the lateral coordinates on the plane are given by 
xii, and Q summarizes the azimuthal (8) as well as the polar {(f)) angle on the sphere. 

Next, the inverse propagator Ai 22 011 the sphere S 2 can be determined from M.22^22 = 
1 which in terms of spherical coordinates reads 



r 2 



[ MoiiC, ft; tt')M 22 {(; ft'; ft") = 4<K0 - $")5(cos6 - cos#"), (13) 



with M.22 given by Eq. ([5]), and J n = J Q 27r d<p /_ x d(cos6>). This equation can be solved for 
M.22 by expanding all quantities in spherical harmonics YJ m (ft). As mentioned above, only 
the monopole order I = = m is relevant here, yielding 

mkq = r^nrrw < 14 > 

monopole 47r r z smh(r | £| ) 

The last propagator in the chain (Eq.([9])) evaluated on the surfaces Si and 5*2 reads, again 
in propertime representation, 



M 2 l((- ft - X\\) = I ^ C~ C~ 4?r2 [ (r COS ^ S ' n 6 ~ Xl )2+(r S ' n ^ S ' n e - x 2) 2 +(g+r cos ef]T ^ 

Jo y/AnT 

With this, all correlators needed for the calculation of the Casimir energy between the 
sphere and the plane to monopole order are at hand. For the Casimir- Polder limit, it 
suffices to consider the lowest-order term of the n sum in Eq. (Q , as discussed below. This 
term reads 



tr (AM 12 M^M 21 ) = ^ / / / / r 4 AA< 12 (C;^||;ft)A^ 2 ~ 2 1 (C)-M 21 (C;ft , ;5 ? ii 

Z7T J( Jx» Jn Jci' 



(16) 



Let us now rescale all dimensionful quantities by the scale set by the sphere-plate distance 
if, i.e., x\\ — y x\\H, ( — ► jj, and consequently also the propertime parameters S — > and 

T — ► jfi. In the limit ^ < 1, the inverse propagator on the sphere ffT4"l) can be expanded 

as 

M£(C) = ±\(l+0(l\C\)). (17) 



47rr 3 V \H 

Even though ( ranges from to 00, the ( integral receives its dominant contributions from 
( = 0(1), such that the order estimate of Eq. f|T7|) is meaningful. The inverse propagator 
M.22 thus becomes independent of the imaginary frequency to lowest order. 

Collecting all the dimensionful factors from the rescaling, the trace Eq. ([TBI) turns out to 
be of order O (jp) and the remaining dimensionless integrals are a pure function of r/H. 

In the limit r -C H, A.Mi 2 (C; ^||> ft) an d -M 2 i(C; ft'; x\\) become independent of Q and Q', 
such that the two solid-angle integrations just contribute a factor of I671" 2 . We then obtain 
for Eq. (1TB1): 



tr(AM 12 M 22 l M 21 ) =T^2 JJ^ AM 12 ((; ^)M 21 ((; £„) + O (j^J 
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with the associated dimensionless propagators 

-i roc 

7T2 Jo 




1 f°° dT _-4ir 2 (l+3f,)f .._ ( C 2 



M ^ = Vii„ T e ' «' ^{^t)- (19) 

Finally, the remaining integrations in the trace expression can be performed straightfor- 
wardly in the following order: First, we perform the integral over ( which is purely GauBian, 
then the integral over xn is done using polar coordinates. At last, the integration over the 
variables of the proper time representation T and S is performed. The result for the 
lowest-order trace term then reads 

tr (AM 12 M 22 1 M 21 ) = T E i^ + (20) 

As higher orders in the n sum involve more propagators between the sphere and the plate 
and consequently further powers of H in the denominator, only the n = 1 term survives 
in the Casimir- Polder limit. Using Eq. (Q, the energy between the sphere and the plane 
consequently yields in the Casimir-Polder limit r H 

which agrees with [SJ M> US EH E] . 

As we have seen, the spatial integrations over the surface of the sphere S2 become trivial 
in the Casimir-Polder limit. Due to this fact, the integrations over the remaining lateral 
coordinates xu on the plate (cf. Eq. ([15]) ) could also have been performed in momentum 
space since the flat plate itself is translationally invariant along these directions. However, 
in the following section we will extend our investigations to surfaces which are uniaxially 
structured along the direction x\. For this purpose, it is expedient to Fourier transform 
Eqs. (I18j) and f|T9l) only with respect to the 2-component p 2 to momentum space. Due to a 
remaining Lorentz invariance in time and x 2 direction, the integrand only depends on the 



combination of momenta q = y p 2 + ( 2 , such that we obtain 
tr (AM 12 M 22 l M 21 ) = T E — J dq J dx, qAM 12 (q; x^M^q; x t ) + O f^j . (22) 

After the execution of the S and T integrals, the propagators are given by modified Bessel 
functions of the second kind, 

1 q 



AM 12 {q;x x ) = - ^— K^q^/l + xj 
n a/1 + x\ 



M 21 (q; Xl ) = —K [q x /l + xj). (23) 
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This representation is suitable for a generalization to a uniaxially corrugated surface, as is 
done in the next section. 

4 Scalar Casimir-Polder potential for uniaxially cor- 
rugated surfaces 

We now extend the above method to uniaxially arbitrarily corrugated surfaces Si. For 
simplicity, we consider deformations along x\ which can be parametrized by a height func- 
tion h(xi) (overhangs could also be included by an appropriate parametrization). The 
four-vector pointing onto the structured surface reads x = (x ,Xi,X2,h(xi)). The center 
of the sphere is again located at x\ = £2 — 0, £3 = H. Here we write H instead of just 
H to point out that it denotes the position of the sphere at a mean distance H from the 
surface. For corrugated surfaces the actual distance H between the surface and the sphere 
is a function of the direction of corrugation H = H(xi) (cf. also Fig. [TJ). As for the planar 
situation considered earlier, it holds that H — H. 

As the inverse correlator on the sphere does not change and the correlators between the 
surfaces are still easy to compute, the main challenge for extending the previous calculations 
to corrugated surfaces is the determination of the inverse propagator M.^ on the structured 
surface Si. For the planar case, translational invariance along the x directions allows for 
diagonalization in momentum space and thus for an explicit solution as given in Eq.( fl0i) . 
By contrast, a structure on the surface breaks translational invariance in the direction of 
corrugation and the functional inverse of M.n is not known analytically in the general case. 

Let us first derive a suitable representation of the problem. The desired quantity A4ii 
is defined by the equation 

f MniCrfrfM^foxtf') = 6(x'-x%, (24) 

JxeSr 

where Ain((; x'; x) is given by the Fourier transform of Eq. (jSJ) with respect to the time 
coordinate. The integration over the surface in this case is defined by 



(25) 

xz=h(xi) 

with the determinant of the induced metric given by g(xi) = 1 + {d\h{xi)) 2 . Thus Eq. ( 12~4"|) 
explicitly reads 

f Vi(^0jV( 11 (C;4;x||)^lri L (C;^||;4) = -^=5(4 -5[f), (26) 

where x» summarizes the flat coordinates x±,X2 along the two-dimensional surface. 

In principle, we could try to solve this equation numerically and then plug in the solution 
into the our formulas for the Casimir energy. However, due to the singularity structure 



xeSi 



Xi 
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of this equation, and also since we rather need the operator product AM.i 2 = M. n M.\ 2 
than M-il alone, we now multiply this equation with a/ ' g{x'{)M.\ 2 {C > \ ) from the 

left and integrate both sides of Eq. (1261) over the lateral coordinates x'/, . One finds that 

f y/^M 11 (C^fx ll )AM 12 (C ] x r ,x'f;,x , i) =M 12 (C;4;4',x^), (27) 

ii 

where AM. X2 explicitly includes the metric factor of the respective structure. Since the 
above matrix elements are still diagonal in p 2 in our uniaxial setup, we take the Fourier 
transform of Eq. (J27I) with respect to the 2-component. Slightly renaming the coordinates, 
we finally get 

/ V g{xi)Mii((,p 2 ; x[; xi)AMi2((,p 2 ; x x ; x", 4') = Mi 2 {C,p 2 ,x' x ;x'{,xl). (28) 

J X\ 

This equation is indeed more suitable for a numerical integration. As a final step, we again 
go over to dimensionless variables by rescaling all dimensionful quantities by the distance 
parameter H J§ Also introducing the combined momentum q = ->/C 2 + p\ results in 

^/gJs!)Mn(q; x; x)AM 12 (q; x) = M 12 (q; x'), (29) 

where we have dropped the coordinate subscript "1", since x = x\ is the only relevant 
direction in this Green's function problem. The dimensionless propagators M.\ 2 = M. 2 \ 
and M.\\ are given by 



Mn(q;x';x) 



Mi 2 {q;x) -- 
where the metric factor reads 



Once Eq. fT29|) is solved for AAi X2 , the solution can be plugged into the Casimir energy 
formula. The scalar Casimir energy between a sphere and a surface which is uniaxially 
corrugated along x thus yields in the limit r -C H 

2 Depending on the relevant geometric parameters, a rescaling with a different parameter such as a 
height amplitude or the wavelength of the corrugation may also be useful. 



^K { q\j (x' - x) 2 + ( h(x') - ~h{xf 2 



l-K Q [q\l{x'Y+(h{i')-l' 



(30) 
(31) 



l+[0.Mx)), //(•'') \jl>(rll). 



(32) 
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where 

poo poo 

a :=2 dq dx s/gfflqAM 12 {q; x)M 2 i(q; x), (34) 

JO J-oo 

is a dimensionless numerical factor that depends on the geometry parameters of the con- 
figuration (measured in units of H). The Casimir-Polder limit is obtained in the limit of 
the sphere radius r being much smaller than any other scale, r <C H, A, A, . . . , where A 
denotes a typical amplitude of the corrugation and A a typical corrugation wavelength. 
The factor a thus is a function of a = a(A/H,X/H, . . . ), but it is independent of r. 

From a technical perspective, the result of Eqs. (I33"j) . (134"]) is very simple. It should 
be stressed that already the first trace term in the initial Casimir-energy formula (Q 
includes nine integrations for the trace: one over the imaginary frequency and four times 
two integrations over the lateral surface coordinates. Due to the trivial dependency of the 
integrand on the lateral coordinates of the sphere in the Casimir-Polder limit, the number 
of integrations was then reduced by four; moreover, the n sum is just replaced by its 
first term in this limit. The emerging translational invariance vertical to the direction of 
corrugation reduces the number of integrations by another two. Thus - instead of nine - we 
are left with three integrations: two of them are directly visible in Eq. (1331) . the third one is 
needed for the construction of A.M12 as a solution of Eq. (|29|) . These simplifications make 
the Casimir-Polder limit accessible to numerical integration for arbitrary height profiles. 

The resulting two integrals in Eq. (134"]) are both convergent, non-oscillatory and gener- 
ically exhibiting a simple one-peak structure. The treatment of the singularity structure 
in the Green's functions equation (1291) . however, requires some care and is treated in the 
Appendix. 

5 Sinusoidal corrugation 

As a first nontrivial example, let us calculate the scalar Casimir-Polder potential for a 
sinusoidal corrugation, see Fig. [TJ The potential for this structure is given by Eqs. (1331 . 
(134"]) . where we use h(x) = A sin(co>x + (f>) as a height function appearing in the propagators 
Eq. (|3T]) as well as in the surface metric factor of Eq. (132]) . The phase <ft is used to modulate 
the relative position of the structure beneath the sphere, since the latter is always fixed 
at x — 0. The mean distance between surface and atom is denoted by H, whereas H 
characterizes the distance of the atom to the surface along the global surface normal. Hence 
H can be viewed as a function of <fi in our conventions, H = H((p) with H(—tt/2) = H + A 
at the sine minimum and H(ir/2) = H — A at the sine maximum, i.e., H = always 
corresponds to atom-wall contact, where the limit r -C H is implicitly understood. 

As the crucial building block for the Casimir-Polder potential, we solve the Green's 
function equation Eq. f]29]) numerically on a 1-D lattice in x direction. This requires to 
invert the propagator M.\\ on the corrugated surface. Even though the singularity of 
this propagator at coincident points is integrable in the continuum, the discretized version 
needs to deal with this singularity explicitly. This is done by introducing a regularization 
parametrized by a short-distance cutoff e, which can be removed after the continuum limit 
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Figure 1: Sphere of radius r at a mean distance H above a sinusoidally corrugated surface 
with amplitude A and wavelength u — In our conventions, we fix the sphere at the 
lateral coordinate x — 0, and effectively shift the structure function h(x) by varying the 
phase cj). In this plot, we have set = —tt/2. The distance parameter H measures the 
sphere- surface distance along the global normal, such that H = corresponds to sphere- 
surface contact for all values of 0. 

has been taken. Details of how this procedure is implemented numerically are given in 
Appendix [A] 

In the following, we display our results for the Casimir-Polder energy always normalized 
with respect to the planar-surface case (for consistency, the normalization factor is also 
determined numerically) . In this manner, the geometry-induced effects are better visible. 
Furthermore, we expect that these results for the scalar case give a qualitative estimate also 
for the electromagnetic case for which the normalizing prefactor has a different distance 
dependence. 

In Fig. [2j we plot -E' S ine/-E'pianar as a function of the horizontal position of the sphere 
above the sinusoidal corrugation between = — ir and = tt for three different mean 
separations H/A = 4, 2, 1.25. The corrugation frequency is chosen to be cuA—1, all units are 
set by the corrugation amplitude A. As expected, above the corrugation minimum, e.g., at 
x = — 7i / 2, the Casimir-Polder potential lies above the planar estimate since the plate bends 
towards the sphere. Analogously, above the corrugation maximum at tt/2, E sine / '-Epianar < 
1. It is clear that the ratio -E S ine/-E P ianar — > 1 as H/A — > oo, since the corrugation cannot be 
resolved anymore for greater distances. However, it is quite noticeable that the deviation 
from the planar case is still up to 10% even at large separations H/A = 4. 

In Fig. El we display E sine /E p i 

anar 

function of the vertical position of the sphere 
above a minimum of the corrugation (0 = — tt/2) for different corrugation frequencies 
ujA = 1, 2, 3. In the limits H/A — > oo and H/A — > 0, we find that -Esine/-Spianar - > 1- This 
is expected, since in the first limit the corrugation of the plate cannot be resolved as it is too 
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^sine/^planar 




X 



Figure 2: We give the ratio E s i ne / E v \ anax for three different separations H/A = 4 (blue) 
H/A = 2 (red) and H/A = 1.25 (green), from bottom to top. The frequency of the sine 
structure is kept fixed at uoA = 1. Above the corrugation trough at <p = — |, it holds 
that i?sino/-Spianar > 1, since the influence of the concavities of the corrugation towards 
the sphere is not accounted for by -Epianar- Similarly, one finds E S i ne /E p i anar < 1 above the 
maximum of the corrugation at = |. For larger ratios of H/A, i.e. larger distances, the 
result approaches the planar case. Note that even at larger separations H/A = 4, a pure 
planar approximation deviates from its true value by up to 10%. 



small compared to the distance. In the second limit, the corrugation is irrelevantly large 
compared to the distance, i.e. the sphere does not notice it locally. In the region where 
H ~ A, the potential for the corrugated surface clearly deviates from the corresponding 
planar case. One can see that the effect becomes more pronounced for larger corrugation 
frequencies, i.e. shorter surface periodicity. 

We identify various regimes which can be classified in terms of an anomalous dimension 
r) which measures the deviation of the Casimir-Polder potential from the planar case, 

-^corrugation ~ JJ^+V ' (^) 

with rj = for the planar case. At small distances, H/A < 1, we find a linear increase 
of the normalized potential E sine /E p \ anar with H/A, implying an anomalous dimension 
of T] = — 1. A linear fit to the short-distance data (not shown in Fig. [3]) in the well 
yields -E'sinc/'-E'pianar ~ 1 + (3{H/A). The linear coefficient (3 depends on the frequency, 
j3 = (3{uj/A) ~ 0.5, 2.3, 5.2 for ujA = 1, 2, 3; within the numerical accuracy, this dependence 
is compatible with a power law (3 ~ (u)A) 2 . 
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At larger distances H/A ~ the normalized energy develops a peak. Various 

regimes can be identified near the peak and also in the drop-off region. The increase 
towards the peak as well as the decrease right beyond the peak can be characterized by 
power laws parametrized by an tu-dependent anomalous dimension. Towards the peak, 
we find r] ~ —0.33, —0.57, —0.67 for ujA = 1, 2, 3, and the fit beyond the peak yields 
i] ~ 0.4, 1.0, 1.6 for ujA = 1, 2, 3. For even larger distances near H/A ~ 10, we observe that 
all normalized energies approach a universal curve being characterized by an anomalous 
dimension rj = 0.2; in particular, the anomalous dimension shows no sizeable uj dependence 
anymore. 

Whereas this observation might come as a surprise in the present formalism, it can 
easily be interpreted in the framework of the worldline picture of the quantum vacuum 
[T2] . In this picture, quantum fluctuations are mapped onto random paths characterizing 
the spacetime trajectories of these fluctuations. In order to contribute to the Casimir 
interaction energy, such a trajectory has to intersect with both surfaces, the sphere and the 
corrugated plate in the present case. This implies that the fluctuation has an average extent 
of the order of the surface separation H. Due to isotropy of the vacuum fluctuations, the 
relevant worldlines also have a lateral extent of this order. This implies that the fluctuation 
integral also averages over structures of the corrugation which are smaller than H. Higher 
corrugation frequencies with ujH 3> 1 thereby become irrelevant for the Casimir-Polder 
potential, as is demonstrated by the universal drop-off for larger H/A. 

For even larger distances H/A ^> 10, the power law cannot continue for arbitrar- 
ily large H/A, since the Casimir-Polder potential eventually has to approach the planar 
limit. In this large-distance regime, we have only a few reliable numerical data points, 
as the discretization artifacts increase, once the lattice spacing approaches the corruga- 
tion wavelength. The available data is compatible with a logarithmic approach towards 

^sine /^planar ~> 1 for H/A -»• OO. 

Finally, we compute the Casimir-Polder potential above a maximum of the sine struc- 
ture at cf) = +7r/2. As expected, the Casimir-Polder energy is always smaller than in the 
planar the surfaces bends away from the atom and approaches the planar result 

in the two limits H/A — > and H/A — > oo, see Fig. HJ Starting from an initial decrease 
of the normalized energy for small distances H/A, a power-law decrease develops towards 
the dip with rj ~ 0.09, 0.11, 0.11 for uiA = 1, 2, 3. Beyond the dip near H/A ~ 1, a power- 
law increase follows with anomalous dimension rj = —0.13, —0.16, —0.19 for uA = 1,2,3, 
respectively. Again, we observe a linear u dependence of rj in this regime. Also, a second 
power-law regime is found for larger distances H/A > 10 with an anomalous dimension 
Tj = —0.07 for the uoA = 1 data. Due to an increase of the discretization artifacts, no 
reliable data for larger uj is available, such that the expected universality in this distance 
regime still needs to be shown. 
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Figure 3: Normalized Casimir-Polder energy -E'sine/-E'pianar above a corrugation minimum 
4> = — 7r/2 versus the normalized distance x = H/A for three different corrugation frequen- 
cies io A = 1,2,3, respectively. All units are set by the corrugation amplitude A. Small 
distances are governed by a linear increase with anomalous dimension rj = — 1, cf. Eq. (1351) . 
The drop-off beyond the peak is characterized by an tu-dependent anomalous dimension 
r) ~ 0.4,1.0,1.6 for ujA = 1,2,3. At larger distances H/A ~ 10, all normalized energies 
approach a universal curve with rj ~ 0.2. The corresponding fit regions are indicated by 
horizontal arrows. Also the increase towards the peak can be parametrized by a power-law 
with anomalous dimensions rj ~ —0.33,-0.57,-0.67 for uoA = 1,2,3 (not shown in the 
plot). 



6 Sawtooth corrugation 

As a second example, we study the Casimir-Polder potential for a sawtooth structure, 
where the wavelength A is 2.8 in terms of the amplitude A, i.e. the dominant frequency 
of its Fourier decomposition is ujA ~ 0.45. These parameters reflect the specifications of 
a sawtooth structure used in a recent experiment |34j . For practical purposes, we actually 
use a smoothed, continuous sawtooth-like structure function with wavelength A, starting 
at h(0) = 0, rising linearly to its maximum amplitude A at h(0.8X) and dropping linearly 
to zero again at h(X) = h(0). 

In Fig. 0, we plot -E saw tooth /-^planar above the corrugation minimum. Qualitatively, the 
result is similar to the sine structure and reveals the various analogous regimes. Quantita- 
tively, the peak and consequently some of the anomalous dimensions are more pronounced. 
The increase towards the peak follows a power-law with anomalous dimension rj ~ —0.3. 
For the decrease right beyond the peak at H/A > 1, we find an anomalous dimension 
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x = H/A 

Figure 4: Normalized Casimir-Polder energy £? S i ne / E p \ axiax above a corrugation maximum 
= 7r/2 versus the normalized distance x = H/A for three different corrugation fre- 
quencies ujA = 1, 2, 3, respectively. All units are set by the corrugation amplitude A. 
The increase beyond the dip is characterized by an tu-dependent anomalous dimension 
f] ~ —0.13, —0.16, —0.19 for ujA = 1, 2, 3. At larger distances H/A > 10, a power law with 
i] ~ —0.07 is observed for the us A = 1 curve. The corresponding fit regions are indicated by 
horizontal arrows. Also the decrease towards the dip can be parametrized by a power-law 
with anomalous dimensions rj ~ 0.09, 0.11, 0.11 for ujA = 1, 2, 3 (not shown in the plot). 

of 7] ~ 1.1. At larger distances H/A ~ 10, we again observe a second power law with 
anomalous dimension rj ~ 0.2 which agrees quantitatively with the anomalous dimension 
in the sinusoidal case. 

Within the worldline picture of quantum fluctuations discussed above, this agreement 
can immediately be understood from the fact that the fluctuation integrals again average 
over the corrugation structures small compared to the distance parameter H. 

7 Conclusions 

In this work, we have presented a new approach to Casimir-Polder forces for corrugated 
surfaces which does not rely on a perturbative ordering of length scales. Our approach is 
based on the constrained-functional-integral approach [22] which - for uniaxial corrugations 
- boils down to a one-dimensional Green's function equation along the corrugation direc- 
tion. This equation is accessible to direct numerical integration provided the integrable 
singularity structures are carefully taken into account. For simplicity, we have studied the 
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Figure 5: Normalized Casimir- Polder energy -E'sawtooth/E'pianar above a corrugation min- 
imum of the saw-tooth profile (see text) versus the normalized distance x = H/A for 
ujA ~ 0.45. All units are set by the corrugation amplitude A. The increase towards the 
peak is governed by a power-law with anomalous dimension rj = —0.3. The drop-off be- 
yond the peak is characterized by an (^-dependent anomalous dimension 77 ~ 1.1. At larger 
distances H/A ~ 10, the normalized energy approaches a curve similar to the curves for 
the sine structure with rj ~ 0.2. The corresponding fit regions are indicated by horizontal 
arrows. 



Dirichlet scalar analogue of the electromagnetic Casimir-Polder case, defined by a fluctu- 
ating scalar field satisfying Dirichlet boundary conditions on the surface and the "atom"; 
the latter is modeled by a small sphere in the limit of zero radius in our approach. 

Two periodic corrugations, a sine and a sawtooth function, are studied in detail, re- 
vealing various regimes with distinct distance dependencies. For a classification of these 
distance dependencies, we have introduced the notion of an anomalous dimension charac- 
terizing the deviation of the distance power law of the Casimir-Polder potential from the 
planar case. In particular, the larger-distance behavior H/A > 1 exhibits two different 
power-law behaviors with positive anomalous dimensions above a corrugation minimum 
and with negative anomalous dimensions above a corrugation maximum. In either case, 
the intermediate distance behavior near H/A ~ 1 is characterized by an anomalous dimen- 
sion, the modulus of which increases linearly with the corrugation frequency. 

Most importantly, we have identified a larger- distance regime near H/A ~ 10 where all 
data above a corrugation minimum is characterized by a universal anomalous dimension 
i] ~ 0.2. This anomalous dimension still depends on the position above the corrugation, 
e.g., rj ~ 0.07 above a minimum, but we have found no dependence neither on the shape 
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of the periodic profile nor on the frequency as long as luH ^> 1. Within the worldline 
picture of the quantum vacuum where fluctuation averages are mapped onto random-path 
averages, this universality can be understood from the fact that small-scale structures are 
averaged out and become irrelevant at larger distances. This observation also justifies to 
use the notion of universality and anomalous dimensions, since the fluctuation averages are 
reminiscent to those of critical phenomena. The resulting Casimir-Polder potential can be 
viewed as a "renormalized" effective Hamiltonian where the running IR cutoff is provided 
by the atom-wall distance. 

We would like to stress that this universality as well as the nontrivial power-law behavior 
cannot be deduced from a perturbative analysis, since perturbation theory in the height 
profile is a Taylor expansion in powers of A/H, whereas a nontrivial anomalous dimension 
t] G R corresponds to a (A/H) 71 dependence of the Casimir-Polder law. Therefore, the 
development and use of a nonperturbative method was absolutely crucial for this work. 

Even though we only considered the Dirichlet scalar case, the constrained-functional- 
integral formalism can straightforwardly be extended to the electromagnetic well 
[28] which carries over to a straightforward generalization of our techniques. Since no 
monopole fluctuations exist in the electromagnetic case due to charge conservation, the 
leading-order Casimir-Polder potential in the planar case follows a ~ 1/H A distance law 
instead of ~ 1/H 2 in the present case. But apart from this, we do not expect further 
dramatic differences as far as the corrugation-dependencies are concerned. Therefore, the 
Dirichlet scalar case may be taken as a rough qualitative estimate also for the electromag- 
netic case; in particular, we expect the occurrence of anomalous dimensions of the same 
order of magnitude. 

The notion of anomalous dimension is also of direct use for Casimir-Polder experiments 
based on quantum reflection such as the atom-beam spin echo technique [321 ESI In 
a certain sense, such an experiment measures the local shape of a potential and thus is 
directly sensitive to anomalous dimensions. Indeed, the results of recent measurements 
with corrugated surfaces can be parametrized by anomalous dimensions of order one. Of 
course, a direct comparison between our results and those of an experiment requires much 
more than the computation of anomalous dimensions, since the atoms near the wall can 
move into all directions and not only along the global normal. The full Casimir-Polder 
potential needs to be mapped out, and the time-dependent quantum reflection problem in 
this potential has to be solved. In any case, the approach presented here lays the foundation 
for this future program. 
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A Numerical procedure 



In the following, we detail our implementation for the numerical evaluation of the Casimir- 
Polder potential for arbitrary uniaxial corrugations, cf. sections |5] and El We proceed as 
follows: First, we solve the Green's function equation for the associated propagator AM.12 
Eq. (l29l) by discretizing the equation with respect to the spatially lateral coordinate x. The 
result is then plugged into Eq. (|33|) . yielding the Casimir- Polder energy upon integration 
of x and q. 

For the first step, we introduce two parameters: ±L X which labels the left and right 
cutoff of the spatial integration, and N x denoting the number of spatial discretization 
sites, respectively. In the end, we remove the discretization by a continuum extrapolation 
N x — ■> 00. 

In principle, L x is a physical parameter encoding the physical size of the surface. Here, 
we will not make use of this option of studying finite-size effects, but compute the Casimir- 
Polder potential in the ideal infinite surface limit by extrapolating to L x — > 00. For this, 
we fix the position of the sphere above the plate at x = and choose a symmetric cutoff for 
x G [—L X ,L X ]. The two limits, continuum (N x — > 00) and infinite-length (L x — > 00) limit, 
have to be taken such that the lattice spacing a x = 2L X /N X also goes to zero, a x — > 0. This 
can be ensured by choosing a suitable function L x = L X (N X ), satisfying L X (N X — > 00) — > 00 
and L X (N X )/N X — > as N x — > 00. In practice, we use 

L x {N x ) = ^y/N x N 0x , (36) 

where a,Q X defines a reference lattice spacing at a reference site number N x = N 0x . Note 
that the lattice spacing a x = a x (N x ) = 2L X (N X ) / N x = cl 0x ^JNq x /N x goes to zero in the 
continuum limit N x — ► 00, while L x — > 00 approaches the infinite length limit. Therefore 
all these idealized limits are controlled by one parameter: N x . In practice, the finite- length 
corrections have always been found to be small compared to discretization effects. In 
general, it suffices to choose the reference lattice spacing such that typically L x (N 0x ) = 2H, 
where Nq x specifies the coarsest lattice in the calculation. 

One serious complication arises when discretizing Eq. (1291) : due to the pole of the zeroth 
Bessel function K at its origin, the matrix M 1 ^ that emerges upon the discretization of 
the spatial arguments diverges in its diagonal entries, i.e., for the case when the spatial 
discretization sites lie on top of each other. Whereas these divergencies are integrable when 
solving the problem in the continuum, the discretized matrix becomes singular. Therefore, 
a regularization procedure is required that facilitates to first take the continuum limit 
before the regulator can safely be removed. Here, we use a UV regularization for the 
propagator in Eq. fl3"Ul) for small arguments z controlled by a small parameter e: 

M n (z) = {^ K0{z) > Z ~\ (37) 

1U J |__L(l n ( 2 + e )_ir (e)-ln(2e)) ,z>e' 1 ' 

where z summarizes all arguments of the propagator including both spatial and momentum 
contributions, entering the Bessel function as a single argument, cf. Eq. ( 130]) . The physical 
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result is expected to arise in the limits N x — > oo and e — > with the continuum limit to 
be taken first before the regulator is removed. 

In a numerical calculation where N x and e are always finite, the order of limits done 
by extrapolation requires a careful choice of N x and e. It is already intuitively clear 
that smaller values of e require larger values of N x , since the proper resolution of a more 
pronounced singularity for smaller e needs a finer lattice. As the pole in the inverse 
propagator on the corrugated surface Si persists irrespectively of the corrugation, the 
numerical discretization and regularization errors can be tested in the planar situation 
where the analytical result is known (cf. section [3]): there, the dimensionless factor a 
amounts to j-. 

In Fig. [6j we plot a as a function of the inverse number of discretization sites 1/N X 
for different values of the cutoff e in the planar case. The values for a depend linearly on 
1/N X to a good approximation and appear to converge for different cutoffs e as 1/N X — > 0. 

a 



0.077 - 



0.076 - 



0.075 - 



0.073 - 




Figure 6: Numerical values for the dimensionless contribution of the lowest-order trace term 
a (I34p in the plane-sphere configuration as a function of the inverse number of discretization 
sites N x for five values of the cutoff parameter e, e = 5- 1(T 5 , 2- 1(T 4 , 4- 1(T 4 , 6- 10~ 4 , 8- 10~ 4 
from bottom to top. The analytical value for a is 1/(4-71") ~ 0.07958. For fixed e, the 
result scales linearly with the discretization 1/N X to a good approximation and appears 
to converge with l/N x — > 0, but it is also visible that the gradients of the curves grow as 
e -> 0. 

Next, we extrapolate the values for a linearly to 1/N X — 0; as the linearity persists to 
a good approximation for all values of N x in Fig. [61 it suffices to use only two data points 
for the extrapolation. We give two separate extrapolations for N x = 80, N x = 100 and 
N x = 180 , N x = 200, respectively. The result is plotted as a function of e in Fig. [3 Recall 
that the analytical value for a yields l/(47r) ~ 0.07958 for the flat plate, which is chosen 
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to be exactly the origin of the coordinate system in Fig. [7J 
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Figure 7: Continuum limit for a as defined in Eq. (I34j) after linear extrapolation to 
1/N X -»• (using 1/N X = 80 and 1/N X = 100 in the lower (red) and 1/N X = 180 and 
1/N X = 200 data in the upper (blue) curve) as a function of the cutoff parameter e for values 
of 10~ 4 < e < 0.1. The intersection between the two plot axes is chosen at a — j-, which is 
the exact value for a in the planar case. The two curves can both be separated into several 
regimes. Consider e.g. the lower curve: It holds that for very small values of e (e J$ 0.0045), 
the result for a lies below the analytical value and appears to diverge as e — > 0, since the 
integrable singularity is not resolved by the number of sites N x used for the extrapolation. 
In the region 0.01 J$ e J$ 0.04, a(e) grows linearly with e (for e <J 0.04, we identify a 
a(e) ~ e 2 dependence). An extrapolation a(e — > 0) in the range 0.01 ~ e ;$ 0.04 thus 
provides an estimate for the cutoff- independent value of a (cf. Eq. f l38|) ). For the upper, 
blue curve the respective regimes can also be identified. However, since the resolution of 
the structure is increased by a greater number of N x , the important linear regime is shifted 
to lower e- values. 



The graphs in Fig. [7J can in fact be divided into several regions. Consider, e.g. the 
lower curve: For values of e ^ 0.0045, the extrapolation l/N x — > underestimates a and 
even appears to diverge as e — > 0. This agrees with our expectation that the integrable 
singularity in the Green's function equation has not been properly resolved with the un- 
derlying discretization; higher values of N x would be required for a more reliable estimate. 
This small-e branch therefore corresponds to a region in parameter space where the result 
arising from the correct order of limits (first N x — > oo, then e — > 0) is not yet visible. 

At about e ^ 0.01, a exhibits a clear linear growth with e. For e ^ 0.04, higher power 
corrections become visible. We conclude that the cutoff-dependent factor a(e) can well be 
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approximated by a power series above the value of e ^ 0.01, 

a(e) = «o + a\t + o^e 2 + . . . . (38) 

Thus, by extrapolating the values for a to e = in the region where a grows linearly with 
e, we obtain a cutoff-independent result «o- As for the extrapolation l/N x — > 0, it suffices 
to use only two sites in e in the linear regime to extract «o; of course, also more data points 
for a higher polynomial fit could easily be employed at the expense of computing time. 

From Fig. El we identify for the 1/N X = 80 and 1/N X = 100 data 0.01 < e < 0.04 as 
the region where a(e) grows linearly with e with only very small higher-power corrections. 
Choosing the data points at e = 0.02 and e = 0.025 for a linear extrapolation, we obtain 
a = 0.07970 which nicely matches the analytical value, the error being below 1%@ As a 
check of the continuum limit, an extrapolation using N x = 180 and N x = 200 sites (upper 
curve in Fig. 0,again at e = 0.02 and e = 0.025) yields a = 0.0799554, which is also within 
1% of the analytical value. The small deviations between these two results can be taken as 
a measure for the overall numerical uncertainty. One can see, that choosing larger values 
of N x for the continuum extrapolation also results in an extension of the linear e regime 
to smaller e values. 

It should be mentioned that the choice of required N x values also depends on the 
corrugation parameters. For instance for high values of the corrugation frequency, a better 
resolution is needed; as a rule of thumb, the lattice spacing a x should always be smaller 
than the smallest dominant wave length of the corrugation. 

All numerical calculations for this work have been performed on a standard desktop 
computer with standard linear algebra packages. Depending on the discretization, the 
calculation of a typical data point including continuum limit and regulator removal takes 
on the order of seconds to several minutes. Since the linear-algebra routines scale with 
~ N x , the computational cost for very fine discretizations can rapidly increase. 
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